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Abstract 

Monte-Carlo simulations of abelian projection of T / pure lattice QCD 
show that 1) Polyakov loops written in terms of abelian link fields alone 
play a role of an order parameter of deconfinement transition, 2) the abelian 
Polyakov loops are decomposed into contributions from Dirac strings of 
monopoles and from photons, 3) vanishing of the abelian Polyakov loops in 
the confinement phase is due to the Dirac strings alone and the photons give 
a finite contribution in both phases. Moreover, these results appear to hold 
good with any abelian projection as seen from the studies in the maximally 
abelian gauge and in various unitary gauges. 
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I. INTRODUCTION 



Color confinement mechanism in QCD is still to be understood. Many works |l -13 



have been done to clarify the confinement mechanism on the basis of the idea of abelian 



projection of QCD [14]. The abelian projection of QCD is to extract an abelian theory 
performing a partial gauge-fixing. After an abelian projection, SU(3) QCD can be regarded 
as a U(l) x U(l) abelian gauge theory with magnetic monopoles and electric charges. 'tHooft 
conjectured that the condensation of the abelian monopoles is the confinement mechanism 
in QCD 0. 

An effective U(l) monopole action is derived from vacuum configurations in SU{2) 
QCD after an abelian projection in a special gauge called maximally abelian (MA) gauge 
PJ4] JT0|JT5|1 . Entropy dominance over energy of the monopole loops, i.e., condensation of 
the monopole loops seems to occur always (for all values of the coupling constant f3) in 
the infinite- volume limit when extended monopoles |16| are considered PJ4],|10|1. The abelian 
charge is confined due to the monopole condensation. The confinement of the abelian charge 
after abelian projection means color confinement as shown in [jHf . Monopole condensation 
is known to be the confinement mechanism also in lattice compact QED [|17| -|22| . 

The string tension is a key quantity of confinement. It vanishes at the deconfinement 
transition temperature T c |2"3]] . It was shown that the same string tension can be derived 
from Wilson loops written in terms of abelian link fields alone after the abelian projection in 
the MA gauge Hljpifl . The string tension derived from the abelian Wilson loops also vanishes 
at T c HH| . Moreover, the abelian Wilson loops can be expressed by a product of monopole 
and photon contributions ||@. The monopoles alone are responsible for the string tension 
in T = and T / SU(2) QCD. The same results are obtained also in SU(3) 
QCD H . 



These results strongly support the 'tHooft conjecture |14| in SU(2) QCD. However, 
the above results are restricted to the special MA gauge. What happens in other abelian 
projections is not known yet. For example, abelian Wilson loops in some unitary gauges are 
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too small to derive reliably the string tension [ 24[] , although it does not mean that the string 
tension can not be derived in the gauges. 

A Polyakov loop is another good order parameter, the vanishing of which means color- 
flux squeezing in the confinement phase. It is the aim of this note to show that 1) a Polyakov 
loop written in terms of abelian link fields alone is a good order parameter , 2) it can be 
written by a product of contributions from Dirac strings of monopoles and from photons, 
3) the former alone vanishes in the confinement phase and 4) these results are independent 
of the gauge choice of abelian projection. The first statement in the MA gauge was already 
shown in 



II. DIRAC STRING AND PHOTON CONTRIBUTIONS TO POLYAKOV LOOPS 

We adopt the usual SU (2) Wilson action. To study gauge dependence, we consider here 
three types of abelian projection, i.e., the MA gauge and two unitary gauges. The MA gauge 
is given [|27P by performing a local gauge transformation V(s) such that 

fl = X)Tr(a-3&(5,/x)(73& t (s,/x)) 

is maximized. Then a matrix 

*i00 = E fj)<r 3 W(3, n) + U\s- A, fi)a 3 U(s - p,, fj)) (1) 
is diagonalized. Here 

U(s,fi) = V(s)U(s,fi)V- 1 (s + fi). (2) 

Two unitary gauges considered here are defined by performing a local gauge transforma- 
tion V(s) such that one of the following two matrices is diagonalized: 

X 2 (s) = ]JU(s + (i- 1)4,4) (3) 
Ms) = E U(s, fi)U(s + fi, u)W(s + z>, fj,)W(s t u), (4) 
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where the sum in is over all plaquette directions. We call the former (the latter) Polyakov 
(F12) gauge. 

After the gauge fixing is over, there still remains a U(l) symmetry. We can extract an 
abelian link gauge variable from the SU(2) ones as follows; 

U(s,/j,) = A(s,n)u(s,fj,), (5) 

where u(s, /i) is a diagonal abelian gauge field and A(s, fi) has off-diagonal components 
corresponding to charged matters. Note that a U(l) invariant quantity written in terms of 
the abelian link variables u(s,/i) after an abelian projection is SU(2) invariant PJ26|. 



Now let us show that an abelian Polyakov loop operator after the abelian projection 

v 4 

P = Re[exp{* £ J 4 (s + (i - l)4)0 4 (s + 1)4)}], (6) 

i=l 

is given by a product of monopole and photon contributions. Here Ja(s) is an external current 
taking +1 along the straight line in the fourth direction and ^(s) is an angle variable defined 
from u(s,n) as follows: 



u(s,fi) 



i e~ i0 ^ j 



(7) 



Using the definition of a plaquette variable /^(s) = d^9 u (s) — d v 6^{s) where is a 
forward difference, we get 

Oa(s) = - £ D( fl - s')[d' v U(s') + d^&Ms'))}, (8) 

s' 

where D(s — s') is the lattice Coulomb propagator and d' u is a backward difference. We have 
used d u d' u D(s — s') = —S ss '. Since d' A J^s) = 0, the second term in the right-hand side of (§) 
does not contribute to the abelian Polyakov loop (H). Hence we get 

N 4 

P = Re[exp{-* ]T J 4 (s + (i - 1)4) £ D(s + (i - 1)4 - s')d'J uA (s')}]. (9) 

i=l s' 

The gauge plaquette variable can be decomposed into two terms: 
Uu(s) = f^u(s) + 2vrn^(s), (-7T < f^s) < tt), 



where d 1 f^s) includes only a photon field and n^s) is an integer- valued plaquette variable 



denoting the number of Dirac strings through the plaquette coming out of monopoles |T9 
Hence we get 

P = Re[P 1 -P 2 ], (10) 

N 4 

P x = cxp{-z £ J A {s + (i - 1)4) D(s + 1)4 - s')d'M(s')}, (11) 

i=l s' 

P 2 = exp{-2vrz £ J 4 (s + (i - 1)4) £ D(s + (i - 1)4 - sO^n^s')}- (12) 

i=l s' 

We observe the photon (P p ) and the Dirac-string (P m ) contributions separately: 

P p = Re[Pi] and P m = Re[P 2 ]. (13) 



III. THE VILLAIN FORM OF QED 

The above separation can be done also in the case of compact QED. We first measure 
P p , P m and P adopting the Villain form [2~B| of the partition function on a 8 4 lattice. Since 



there are natural monopoles and DeGrand-Toussaint monopoles in the Villain case of QED, 
we observe P m in terms of the two types of the monopoles. Since the auto-correlation time 
is long for (3 > (3 C , we have to perform Monte-Carlo simulations carefully. We follow the 
same method as done in |15| . The results are shown in Fig. [l|. We find the following. 



1. The monopole Dirac-string data vanish in the confinement phase, whereas the photon 
data remain finite and change gradually for all (3. The characteristic features of the 
Polyakov loops are then due to the behaviors of the Dirac-string contributions alone. 

2. Monopole Polyakov loops show more enhancement than the total ones for (3 > j3 c . 

3. Both types of monopoles give almost the same results. 
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IV. THE MA GAUGE IN SU(2) AND SU(3) QCD 

The Monte-Carlo simulations were done in SU(2) on 16 3 x 4 lattice from (3 = 2.1 to 
j3 = 2.5 in the MA gauge and in the unitary gauges. In SU(3) QCD, we adopted 10 3 x 2 
lattice from f3 = 5.07 to (3 = 5.12. All measurements were done every 50 sweeps (40 sweeps 
in the SU(3) case) after a thermalization of 2000 sweeps. We took 50 configurations totally 
for measurements. The gauge-fixing criterion in the MA gauge is the same as done in Ref. 
2|. 

The results in the MA gauge are shown in the following: 

1. We plot the SU (2) data in the MA gauge in Fig. The abelian Polyakov loops remain 
zero in the confinement phase, whereas they begin to rise from the critical temperature 
(3 C = 2.298 [|30| This was observed already in ||26|1 . It is interesting that the Dirac-string 



contribution shows similar behaviors more drastically. It is zero for (3 < (3 C , whereas it 
begins to rise rapidly and it reaches ~ 1.0 for large (3. On the other hand, the photon 
part has a finite contribution for both phases and it changes only slightly. Character- 
istic behaviors of the abelian Polyakov loops as an order parameter of deconfinement 
transition are then explained by the Dirac-string part of monopoles alone. This is 
consistent with the results in P,|9|,p4]] stating monopoles alone are responsible for the 
value of the string tension. 

2. The same results are obtained also in pure SU(3) QCD in the MA gauge as shown in 
Fig. |||. The monopole Dirac string alone is seen to be responsible for the flux squeezing. 
There is a clear hysteresis behavior showing the first order transition. 

V. THE UNITARY GAUGES 

We next study the case in the unitary gauges. 

1. The data in the unitary gauges are plotted in Fig. [| and Fig. [5| It is very interesting 
to see that qualitative features are similar to those in the MA gauge. Namely, the 



abelian and the Dirac-string Polyakov loops are zero in the confinement phase, which 
suggests occurrence of flux squeezing in the unitary gauges, too. They show finite 
contribution above the critical temperature f3 c . Photon contributions are finite and 
change gradually in both phases. The Dirac-string part of monopoles is responsible for 
the essential features of an order parameter. These are the first phenomena suggesting 
gauge independence of the 'tHooft conjecture. 

2. Comparing the figures in the unitary gauges, one can see both are almost equal, al- 
though the gauge fixing conditions are quite different. The finite values in the decon- 
finement phase in Fig. f| and Fig. [| are much smaller than those in Fig. Q in the MA 



gauge. It is similar to the behaviors of Wilson loops in these gauges studied in |24 
Namely, abelian Wilson loops in the MA gauge enhance drastically and then we could 
determine the string tension. On the other hand, the Wilson loops in the unitary 
gauges are too small to fix the string tension reliably. 

3. There is a problem in extracting monopoles in the unitary gauges in comparison with 
the case in the MA gauge, as seen from the histogram of f^ in some configurations. 
The data in both gauges are shown in Fig. ^| and Fig. |7[ The data in the MA gauge was 
first studied in [31]. In the MA gauge, quantum fluctuation is small and it is reliable 



to separate Dirac strings out of f^ u . However, in the unitary gauges, the separation 
between Dirac strings and quantum fluctuation is more ambiguous as seen from the 
bump around ±7r in the histogram. To be noted, the bump becomes smaller as (5 
becomes larger. This is because quantum fluctuation becomes smaller for larger (3. 
Hence, we may expect that the correct string tension can be derived even in unitary 
gauges when we consider the T = case for large (5 on larger lattices. 
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VI. CONCLUSION AND REMARKS 



In conclusion, our analyses done here strongly suggest that abelian monopoles are respon- 
sible for confinement in SU (2) QCD and condensation of the monopoles is the confinement 
mechanism. We have found first the data suggesting gauge independence of the 'tHooft 
conjecture. 

Our data show the following picture of color confinement due to monopole condensation. 
Choose any U(l) out of SU(2) through an abelian projection. Then quarks and gluons 
behave like charged matters with respect to the U(l) symmetry after the abelian projec- 
tion. There are always monopoles with magnetic charges with respect to the magnetic U(l) 
symmetry dual to the U(l) chosen. The electric U(l) charge is confined due to the dual 
Meissner effect caused by the condensation of the corresponding monopoles. The electric 
charge confinement after abelian projection is equivalent to color confinement as proved in 



14}| . Gauge independence of the confinement mechanism appears in this way. 

In p4 |, abelian Wilson loops have been measured in various gauges. The abelian Wilson 
loops in the unitary gauges do not show abelian dominance and take a similar value to 
that without gauge-fixing. Suggested from the study, we have tried to measure abelian 
Polyakov loops without gauge-fixing, although such quantities are gauge variant. Abelian 
link fields are defined by choosing any one of isospin directions. Surprisingly enough, we 
have obtained a similar behavior as shown in Fig. |8|. The abelian Polyakov loop is zero 
in the confinement phase and shows rising at the critical (3 C . The responsibility of the 
Dirac string is also seen. The abelian quantity without gauge-fixing is variant under local 
SU(2) transformation. Hence to see whether the above behavior includes some physical 
meaning or not, we have also measured abelian Polyakov loops composed of abelian links 
defined randomly at each site. Namely abelian links in the time direction at each site are 
defined to take a random isospin direction. The data are shown in Fig. |9|. Qualitatively 
similar behaviors are obtained, although the finite values of the total and the Dirac-string 
contributions in the deconfinement region are smaller. These data may support the above 
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picture of gauge independence of the 'tHooft idea. To prove gauge independence definitely 
is to be studied in future. 

This work is financially supported by JSPS Grant-in Aid for Scientific Research 
(B)(No.06452028). 
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FIG. 1. Monopole Dirac string and photon contributions to Polyakov loops in the Villain 
model of compact QED. 
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FIG. 2. Monopole Dirac string and photon contributions to Polyakov loops in the MA gauge 
in SU(2) QCD. 
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FIG. 3. Monopole Dirac string and photon contributions to Polyakov loops in the MA gauge 
in SU(3) QCD. 
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FIG. 4. Monopole Dirac string and photon contributions to Polyakov loops in the Polyakov 
gauge. 
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FIG. 5. Monopole Dirac string and photon contributions to Polyakov loops in the F12 gauge. 
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FIG. 6. Histogram of in the MA gauge for some f3. 



17 



10000 



SU2: POLYAKOV GAUGE 



beta=2.1 
beta=2.5 



5000- 



°-4 




FIG. 7. Histogram of / M „ in the Polyakov gauge for some (3. 
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FIG. 8. Monopole Dirac string and photon contributions to abelian Polyakov loops without 
gauge- fixing. 
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FIG. 9. Monopole Dirac string and photon contributions to abelian Polyakov loops without 
gauge-fixing. Abelian links are defined randomly at each site. 
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